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Hot subdwarf B stars (sdBs) are extreme horizontal branch stars beheved 
to originate from close binary evolution. Indeed about half of the known sdB 
stars are found in close binaries with periods ranging from a few hours to a 
few days. The enormous mass loss required to remove the hydrogen envelope of 
the red-giant progenitor almost entirely can be explained by common envelope 
ejection. A rare subclass of these binaries are the eclipsing HW Vir binaries 
where the sdB is orbited by a dwarf M star. Here we report the discovery of an 
HW Vir system in the course of the MUCHFUSS project. A most likely substellar 
object (~ 0.068 Mq) was found to orbit the hot subdwarf J08205+0008 with a 
period of 0.096 days. Since the eclipses are total, the system parameters are 
very well constrained. J08205+0008 has the lowest unambiguously measured 
companion mass yet found in a subdwarf B binary. This implies that the most 
likely substellar companion has not only survived the engulfment by the red-giant 
envelope, but also triggered its ejection and enabled the sdB star to form. The 
system provides evidence that brown dwarfs may indeed be able to significantly 
affect late stellar evolution. 

Subject headings: binaries: spectroscopic — stars: individual(SDSS J082053. 53-1-000843. 4) 
— brown dwarfs — stars: horizontal-branch 



1. Introduction 



Hot subdwarf B stars (sdBs) are core helium-burning stars with hydro gen envelope s 
too thin to sustain hydrogen shell burning and have masses of about 0.47 Mq ( lHeberll2009l ). 
The large fraction of close binaries - about half of the known sdB stars are members of 



short-period (P < 10 days) binaries (IMaxted et al.ll200ll : iNapiwotzki et al.ll2004al ) - can be 
explained by binary evolution models. The required extraordinarily large mass loss in the 
red giant phase is triggered by the for mation of a common e nvelope, which is finally ejected. 



Binary population synthesis models ( iHan et al.l |2002| . |2003[ ) are successful in matching the 



observed properties of known systems qualitatively. The existence of apparently single sdB 
stars poses another problem. However, even in this case binary evolution com es to the rescue 
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The existence of eclipsing sdB-|-dM binaries of HWVir type with very short orbital 
periods (0.1 — 0.26 d) and very low companion masses between 0.1 and 0.2 (e.g. 
For et al.ll2010l : l0stensen et al.ll2010[ ) shows that stars close to the nuclear burning limit of 
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0.08 Mq are able to eject a common envelope and form an sdB. Claims have 



3een made that 



the p rimaries of AA Dor flRauchll2000l : iRucinskil l2009f ) and HS 2231+2441 fl0stensen et al. 



20081 ) are low mass sdBs and the companions therefore substellar. However, in the former 
case this has been refu ted by the measurement of the companion's radial vel ocity (RV) curv e 
( IVuckovic et al.l 120081 ) and in the latter case by the gravity measurement (IFor et al.ll2010l ). 
Substellar co mpanions to sdB stars have been found using the light travel time technique 
( ISchuhl l2010l . and references therein). However, these systems have wide orbits and are 
therefore unlikely to have experienced a common envelope phase. None of these companions 
influenced the evolution of its host star, but the presence of such objects in wide orbits may 
be an indication of current or former close substellar companions to sdBs. 

Here we report the discovery of the short-period eclipsing HW Vir type binary J08205+0008 
from the MUCHFUSS project. This system is regarded as the first one in which a close sub- 
stellar companion to an sdB star has been detected unambiguously. 



2. Target selection and radial velocity curve 



The project Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS 
(MUCHFUSS) aims at finding sdBs with compact companions like massive white dwarfs 
(WDs, M > 1.0 Mq), neutr on stars or b l ack hol es. Details on the survey and target selection 
procedure are prov ided in iGeier et al.l (j2010al ). and an analysis of seven sdB binaries in 
Geier et al.l ( l2010bl ). The same selection criteria that we applied to find such binaries are 
also well suited to single out hot subdwarf stars with constant high radial velocities in the 
Galactic halo and sea rch for hyperyelocity stars. First results of the search for hypervelocity 
stars are presented in lTillich et al.l (120101 ). 



The MUCHFUSS target selection strategy is tailored to single out RV variations on 
time scales of half an hour or less. Such variations may indicate the presence of short-period 
systems of relatively low RV amplitude or longer-period binaries with high RV amplitudes. 
The latter are the prime targets for the core programme of MUCHFUSS. Obviously, the 
campaign is also bound to find short-period, low RV amplitude systems with low mass 
stellar or even substellar companions. 

SDSS J082053.53+000843.4 (GALEX J082053.6+000843, in short J08205+0008, g = 
14.9 mag ) was classified as an sdB star by colour selection and visual inspection of SDSS 
spectra (lAbazajian et al.ll2009[ ). which are flux calibrated and cover the wavelength range 
from 3800 A to 9200 A with a resolution of R = 1800. The six individual sub-spectra showed 
significant RV variability and the star became a high-priority target for the MUCHFUSS 
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spectroscopic follow-up. Eighteen spectra were taken with the EF0SC2 spectrograph [R ~ 
2200, A = 4450-5110 A) mounted at the ESONTT. Five additional spectra were taken with 
the Goodman spectrograph mounted at the SOAR telescope {R ~ 2500, A = 3500 — 6160 A). 

The radial velocities were measured by fitting a set of mathematica l functions to the 



hydrog en Balmer lines as well as He I lines using the FITSB2 routine (INapiwotzki et al. 



2004b|). The continuum is matched by a polynomial, line wings and line core by a Lorentzian 



and a Gaussian, respectively. The orbital solution (see Table H]) was derived based on th e 
spectra from SDSS and the ones taken with EF0SC2 as described in iGeier et al.l ( l2010bl ). 
The phase coverage of the RV curve is very good (see Fig. [1]). The best fit orbital period is 
0.096 ±0.001 d. Two ahas periods are possible at P=0.088 d and 0.108 d, but can be rejected 
by the analysis of the light curve. 



3. The light curve 



Photometric light curves were obtained on three separate nights between November 
29, 2009 and January 13, 2010, with the Flemish 1.2-m Mercator Telescope on La Palma, 
Canary Islands. During the period that these observations were made, we were in the 
process of commissioning an upgrade to the Merope CCD c amera, fitting it with a large E2V 
frame transfer CCD with 2048x3074 pixels (j0stensenll2OlOl ). Although the upgraded camera 
suffered from some technical issues, and the observing conditions were far from perfect, the 
deep primary and shallow secondary eclipses superimposed on a strong reflection effect are 
perfectly clear in all the light curves (see Fig. |2]). The data were obtained using the Geneva 
R-band filter and exposure times from 16 to 90 seconds. 

The light curve was analysed with the M ORO code, which is based on the Wilson- 
Devinney approach ( iWilson fc Devinneylll97l[ ) but takes into accou nt radiative interaction 
between the components of hot, close binaries (IDrechsel et al.lll995l ). 



The light curve was binned over narrow time intervals. The photometric period was 
determined by measuring the timespans between three consecutive primary eclipses. The 
result 0.096 ±0.001 d is perfectly consistent with the best spectroscopic solution (see Sect. |2]). 
We used Wilson-Devinney mode 2, which poses no restrictions to the system configuration 
and links the luminosity and the temperature of the second component by means of the 
Planck law. The gravity darkening exponents (^'1,2) and t he linear limb darke ning coefficient 
of the sdB primary {xi{R)) were fixed at literature values (IDrechsel et al.ll200ll . and references 
therein), whereas the linear limb darkening coefficient of the secondary {x2{R)) was fixed 
to 1.0 because all light curve solutions converged at this value. The bolometric albedo of 
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the primary was also fixed to 1.0. The temperature of the sdB was taken from the spectral 
analysis (Teff(l) = 26 000 K). 

The remaining adjustable parameters are the inclination, the temperature of the second 
component (Teff(2)), the Roche potentials (^^1,2), the bolometric albedo of the secondary 
{A2), the radiation pressure parameter {61) and the luminosity of the hot component. The 
fractional Roche radii (ri_2) in units of the orbital separation a were calculated using the 
Roche potentials and the mass ratio q. We used the binary mass function derived from 
spectroscopy to calculate possible mass ratios for a range of primary masses. A grid of light 
curve solutions with different mass ratios was calculated. In order to derive errors we created 
500 new datasets with a bootstrapping algorithm by random sampling with replacement from 
the original dataset. In each case a light curve solution was calculated in the way described 
above. The standard deviations of these results were adopted as the error estimates for the 
parameters. 

The flat-bottomed eclipses, and the way the secondary eclipses reach down to the flux 
level just before and after primary eclipse, indicate that the secondary is totally eclipsed by 
the primary. The eclipses ensure that the inclination of the system and the relative radii 
of the components are very well constrained, but the usual degeneracy still remains in the 
mass ratio and the ratio of the components' radii and the size of the orbit. Table [1] shows 
light curve solutions for the most likely sdB masses after combining the photometric and 
spectroscopic analyses. Fig. [2] shows an example model flt to the light curve. 



4. Atmospheric and stellar parameters of the sdB star 



Atmospheric parameters have been determined by fltting a grid of synthetic spectra, cal- 
culated from line-blanketed, solar-metalicity LTE model atmospheres (iHeber et al.ll2000l ). to 
th e hydrogen Balmer a nd helium lines of the SDSS and SOAR spectra in the way described 



m 



Geier et al.l (l2010al ). The single spectra have been corrected for their orbital motion 



and coadded. In order to investigate systematic effects introduced by the individual instru- 
ments, especially the different resolutions and wavelength coverages, the parameters have 
been derived separately from spectra taken with SDSS {S/N = 83, T^s = 26000 ± 1000 K, 
log^ = 5.37 ± 0.14) and SOAR {S/N = 61, T^s = 26900 ± 300 K, log^ = 5.51 ± 0.04), 
respectively. The weighted mean values have been calculated and adopted as flnal solutions 
(see Table [1]). 



The contribution of light from the irradiated surface of the cool companio n in HWVir 
type binaries can lead to systematic shifts in the atmospheric parameters (e.g. iHeber et al. 
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2004; 


For et al. 


2010; 


Miiller et al. 


2010) 



cover most of the orbital phase is not good en ough to resolve t his eff ect. The statistical errors 



are higher than the expected modulations. iDrechsel et al.l ( 120011 ) adopted errors of 900 K 



in Tcff and 0.1 dex in log (7 to account for this effect in the case of the HWVir type binary 
HS 0705+6700. Since the orbital period as well as the atmospheric parameters of this binary 
are similar to the ones of J08205+0008, we adopt similar uncertainties here. Systematic 



errors introduced by different model grids are typically smaller than that (ILisker et al.ll2005 
Geier et al.ll2007l ). 



While the inclination of the system is well constrained by the light curve analysis, 
there remains a degeneracy between the masses and the radii of the components. The 
surface gravity of the sdB determined in the quantitative spectroscopic analysis provides an 
additional constraint since it only depends on the mass and the radius of the subdwarf. To 
proceed, however, we need to constrain the sdB mass from evolutionary models. 

In Fig. [3] we compare the position of the star in the {Tes, log5')-plane to other HW Vir 
stars as well as to two sets of models. The first one represents the canonical picture of EHB 
evolution, while the second one recalls post-RGB evolution, which means that the sdB star 
has left the RGB early and did not ignite helium in the core at all. 

According to Fig. [3] the star is situated on the Extreme Horizontal Branch (EHB) 
consistent with being a core helium-burning star as are the other HW Vir stars (except 
AA Dor). Since the orbital period is short, i t was formed via common envelope ejection. 
Population synthesis models (IHan et al.ll2002l . 120031 ) predict a mass range of M^dB = 0.37 — 
0.48 with a sharp peak at 0.47 M0 for sdBs in binaries formed in this way. Even lower 
sdB masses (down to 0.3 Mq) are possible, when a more massive progenitor star (2 — 3 M©) 
ignites core helium-burning under non-degenerate conditions. Since this formation channel is 
predicted to be rare, the results of the quantitative spectroscopic analysis are fully consistent 
with an EHB model of the canonical m ass. Hot subdwarf masses derived from asteroseismic 
analyses (e.g. Ivan Grootel et al.ll2010[ ) or from analyses of eclipsing binaries are in general 



agreement with this picture. 

In spite of the consistency discussed above, it may still be premature to adopt the 
canonical mass. It has been pointed out that the sdB stars in AA Dor, HS 2333+39 27 and 
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evolution is interrupted by the ejection of a common envelope before the core has reached 
the mass required for helium ignition. 



Driebe et al.l ( 1l998l ) calculated evolutionary tracks of these so-called post-RGB objects. 
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In Fig. [3]we see that the evolutionary track for a 0.25 Mq post-RGB star is also consistent 
with the atmospheric parameters of J08205+0008. The time it takes for such a low mass 
star to cross the region of the Teg — log g plane where canonical subdwarfs spend their time 
is only ^ 2Myr (for the model plotted in Fig. [3]), compared to 100 — 150 Myr for canonical 
EHB stars. Hence, such a low mass solution is unlikely, but cannot be completely ruled out. 
In Tabled] we show solutions for both the canonical EHB and the post-RGB scenario. 

To further constrain the mass of the sdB star we can make use of the spectroscopic 
gravity determination. However, the \ogg derived from the spectra is fully consistent with 
both solutions within the uncertainties. To tighten this constraint, the gravity would need 
to be determined to better precision (Alogf/ ~ 0.05). 



5. Constraining the nature of the companion 



We now turn to the cool companion. Fig. H] shows its mass-radius relation derived 
from the light curve analysi s. Theoreti cal mass-radius relations for low mass stellar and 



substellar objects taken from lBaraffe et al.l ( 120031 ) are given for comparison. Since the normal 
progenitors of sdBs are expected to be stars around 1.0 M©, it seems reasonable to assume 
an age of a few Gyr for the system. The relations intersect at two very different mass ranges 
- a low mass one close to 0.25 Mq and a high mass one at 0.78 Mq. The high mass solution is 
very unlikely, because su ch high sdB masses are neither predicted by theor y (IHan et al.ll2002 



2003 : IZhang et al.l l2009l ) nor have ever been measured empirically (e.g. Ivan Grootel et al. 
2010h . 



For the canonical sdB mass of 0.47 Mq and the corresponding companion mass of 
0.068 Mq, the radius of the companion would have to be ~ 20% larger than expected from 
theory (see Fig. H]). The irradiation by the sdB primary should influence the radius of the 
companion, which may be inflate d. This effect has been measured in the case of hot Jupiter 



(iLittlefair et a 



Baraffe et al. 



planets (e.g. lUdalski et al.ll2008h . the accreting WD+BD binary SPSS 103533.03+055158 4 



20061 ) and the eclipsing MS+BD binary CoRoT-15b f iBouchv et al.l 120111 ). 



2OO3I ) calculated mass-radius relations for planets including the influence of 



irradiation. The effect is of the order of up to 10% for hot Jupiters orbiting solar type stars. 
More massive objects should be less affected. On the other hand the sdB is more luminous 
than a solar type star and the separation of the system is less than one solar radius. Re- 
cently iParsons et al.l (120101 ) modelled the inflation of a late M star irradiated by a hot WD 
(Teff = 57 000 K) in a 0.13 day orbit and showed that the measured 10% increase in radius 
compared to theoretical models can be reproduced in this way. For a discussion of othe r 
possible solutions to this issue (e.g. cold spots on the BD surface) see iBouchy et al.l ( 1201 ll ). 
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For a low sdB mass of ~ 0.25 Mq the companion's mass-radius relation is also consistent 
with theory. In this case the companion mass would be 0.045 M©. A third possibility may 
be that the progenitor of the sdB was originally more massive. In this case the sdB mass 
could be as low as ~ 0.3 Mq and the system as well as the substellar companion would be 
much younger. Since young BDs are considerably larger, this may also lead to a consistent 
solution (see Fig. H]). 

In the cases discussed above, the mass of the companion ranges between 0.045lg og2 ^0 
and 0.068l°:gi M(7). The most con servative theoretical lower limit for core hydrogen-burning 
(~ O.OTMy. lChabrier et al.ll2000l ) is right at the border of this range. We therefore conclude 
that the companion is most likely a brown dwarf. However, given that the sdB mass is not 
strictly constrained, the companion may also be a star of extremely low mass. 



6. Discussion 

We have presented a spectroscopic and photometric analysis of the HW Vir type eclips- 
ing sdB star J08205+0008, discovered in the course of the MUCHFUSS project. The com- 
panion turns out to be a very low-mass, most likely substellar object. 

Although the mass of the sdB is not yet tightly constrained, it is important to stress 
that the companion remains below the core hydrogen-burning limit for reasonable sub dwarf 
masses ranging from 0.25 to 0.47 Mq. The inclination constraint from eclipses means 
that J08205+0008 has the lowest unambiguously measured companion mass yet found in a 
sub dwarf B binary. 

The question whether the sdB is burning helium in its core or not remains open for now. 
Time resolved high resolution spectroscopy is necessary to measure both the f rot sin i and 
the log (7 of the subdwarf with high accuracy. Combined with a high-quality multi-colour 
light curve much tighter constraints could be put on this unique binary system. The fact 
that the sdB is situated on the EHB is a strong argument in favour of the EHB-scenario, 
because post-RGB objects are very rare and not related to the EHB. 

As witnessed by the HWVir type systems, stellar companions with masses as low as 
0.1 Mq are able to eject a common envelope and form an sdB star without being destroyed. 
The case of J08205+0008 demonstrates that even lower mass objects, i.e. s ubstellar ob 



jects, are sufficient. This finding can be used to constrain theoretical models (ISokerlll998 



Nelemans &: Taurid Il998l ) and learn more about the role of substellar companions for the 



formation of single and close binary sdBs. 
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The double-lined spectroscopic WD+BD system WD 0137-349 flMaxted et al.l 120061 ) 
is a binary very similar to J0820+0008, but in a later stage of evolution. It consists of a 
He-core white dwarf of 0.39 orbited by a 0.053 Mq brown dwarf in 0.0803 days. When 
evolving on the white dwarf cooling sequence J08205-I-0008 will therefore appear as a twin 
to WD 0137—349 once it is cooled down to the effective temperature of the latter (15 000 K). 



Based on observations at the La Silla Observatory of the European Southern Obser- 
vatory for programmes number 082.D-0649 and 084.D-0348 and on observations with the 
Southern Astrophysical Research (SOAR) telescope operated by the U.S. National Opti- 
cal Astronomy Observatory (NOAO), the Ministrio da Ciencia e Tecnologia of the Federal 
Republic of Brazil (MOT), the University of North Carolina at Chapel Hill (UNC), and 
Michigan State University (MSU). Based on observations collected with the Flemish 1.2- 
m Mercator Telescope at the Roque de los Muchachos, La Palma, Spain. A.T. and S.G. 
are supported by the Deutsche Forschungsgemeinschaft (DFG) through grants HE1356/45-1 
and HE1356/49-1. R.0. acknowledges funding from the European Research Council under 
the European Community's Seventh Framework Programme (FP7/2007-2013)/ERC grant 
agreement N- 227224 (prosperity), as well as from the Research Council of K.U.Leuven 
grant agreement GOA/2008/04. 
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Fig. 1. — Radial velocity plotted against orbital phase of J08205+0008. The RV data were 
phase folded with the most hkely orbital period. The residuals are plotted below. The 
RVs were measured from spectra obtained with SDSS (rectangles) and ES0-NTT/EF0SC2 
(circles). The errors are formal la uncertainties. 
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Fig. 2. — Phased R-band light curve of J08205+0008. The shght asymetry in the secondary 
echpse is most hkely an artifact caused by the combination of the different datasets. The 
model hght curve for the post-RGB case is overplotted as solid hne. Fits of similar quality 
can be obtained for different mass ratios. 
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Fig. 3. — Tefj— log gf-diagram. The helium main sequence (HeMS) and the EHB band (limited 
by the zero-age EHB, ZAEH B, and the terni i nal-ag e EHB, TAEHB) are superi mposed with 
EHB e volutionary tracks from Dorman et al. ( 1993 ) and a post-RGB track from lPriebe et aP 
f|l998[ ). The position of J08205-I-0008 is indicat ed with a solid diamond. Open diamonds 



mark the posit i on of other HW Vir-like systems ( Charpinet et al. 
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Fig. 4. — The radius of the cool companion plotted against the mass. The black symbols 
mark the relation derived from the spectroscopic and photometric analysis. The three curves 
mark theoretical relations for substellar objects w ith of 1 (dashed line), 5 (dotted line) and 
10 Gyr (solid line) taken from lBaraffe et al.l ( 120031 ) . The dashed vertical lines mark different 
values of the corresponding sdB masses. 
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Table 1: Parameters of J08205+0008 



SPECTROSCOPIC PARAMETERS 



To 


[HJD] 


2455147.8564 ± 0.0006 


P 


[d] 


0.096 ± 0.001 


7 


[kms~-^] 


9.5 ±1.3 


Ki 


[kms""*"] 


47.4 ±1.9 


f{M) 




0.0011 ±0.0001 




[K] 


26700 ± 1000 


\ogg 




5.48 ±0.10 


logy 




-2.0 ±0.07 


^2,M,dB=0.47 


[Mo] 


0.068^°;:^ 


^2,M,dB=0.25 


[Mo] 


0.0451°:°°^ 



LIGHT CURVE SOLUTION 



reff(l) 

91 

92 

xi{R) 

X2{R) 



[K] 



1.0 
26700 

1.0 
0.32 
0.18 

1.0 







M,dB = O.25M0 


MsdB = O.47M0 


g(=M2/Mi) 




0.181 


0.1438 


i 


[°] 


85.87 ±0.16 


85.83 ±0.19 


q{^M2/M^) 




0.181 


0.1438 




[K] 


2958 ± 207 


2484 ± 230 






1.11 ±0.05 


1.09 ± 0.04 


Qi 




3.687 ±0.026 


3.621 ±0.027 






2.732 ±0.009 


2.470 ±0.007 






0.99983 ± 0.00008 


0.99995 ± 0.00004 






0.0282 ±0.002 


0.030 ±0.002 


ri(pole) 


[a] 


0.277 ±0.002 


0.278 ± 0.002 


ri (point) 


[a] 


0.283 ± 0.002 


0.284 ± 0.002 


ri (side) 


[a] 


0.280 ± 0.002 


0.282 ± 0.002 


ri(back) 


[a] 


0.282 ± 0.002 


0.283 ± 0.002 


r2(pole) 


[a] 


0.136 ±0.001 


0.137 ±0.001 


r2 (point) 


[a] 


0.137 ±0.001 


0.138 ±0.001 


r2(sidc) 


[a] 


0.137 ±0.001 


0.138 ±0.001 


r2(back) 


[a] 


0.140 ±0.001 


0.142 ±0.001 



